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Abstract 

Exact Tolman solutions are used to analyse the implications if the 
galactic number count has a fractal form out to a distance of about 
150 Mpc in a universe which is homogeneous on the large scale. It is 
concluded that such a model requires either a non-linear Hubble law or 
a very low density if galaxies trace the total matter distribution. 

1 Introduction 

This paper is based on joint work with Roy ]i.fcartens and Neil Humphreys. 
More details ape given in three joint paperstE citeDR. There is substan¬ 
tial evidence □□ to support a claim that the number of galaxies within a 
volume defined by a distance y has a fractal distribution, 

iV(y) = Ay-' (1) 

where 0 < ;/ < 3 is the fractal index and A is a constant. How large y 
can be is a matter of dispute but there is some agreement that for y < 
150 Mpc the formula (|^ holds with i/ between 1.5 and 2. 

Unfortunately at this stage in the history of cosmology there is no 
convincing theoretical model which predicts such a distribution. Some 
work has been done recentlyu but it is only a beginning. This is a problem 
because it is notoriously difficult to derive a convincing theoretical model 
from data alone. Statistical procedures are excellent if we know the 
model and wish to determine the parameters. However they become 
suspect when we have to derive the model and the parameters from the 
same data. Accumulation of data and improved statistics u help but 
they do not remove unease and doubt, e.g. is the distribution really 
multi-fractal or is the data set incomplete? We can expect the debate 
to continue. 

Our contribution has been to provide a different approach to the 
problem. We assume; (a) a fractal form for the number counts holds 
to a distance (observer area distance) 150 Mpc from us, (b) space-time 
is spherically symmetric about us, (c) Einstein’s field equations hold, 
(d) the matter distribution can be modelled as dust on a suitable av¬ 
eraging scale, (e) beyond 150 Mpc the universe can be modelled by a 
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Friedmann-Lemaitre-Robertson-Walker (FLRW) geometry, (f) the solu¬ 
tions are regular, and there is no shell crossing or surface density layers 
and (g) the Darmois conditions hold at matching surfaces. A conse¬ 
quence of these assumptions is that we can use a Tolman model within 
150 Mpc. Consistent with the picture we take seriously that the obser¬ 
vations are on the observer’s past light pane. This is often not doae by 
astronomers as pointed out by Ellis eit alliil and Laix and Starkmanu and 
theorists do not always do it either Q, but it has signihcant effects, for 
instance on the measuremeaits of the density power spectrum as pointed 
out by Laix and Starkmanu. Even in the spatially homogeneous FLRW 
models it has apwpi trivial effect leading to an in-homogeneous number 
count formula QlijEJ. 

Observational coordinates liiJ provide a powerful tool for analysis on 
the past lighticone of an observer and hence for discussions of real ob- 
servationsEJLJ. They are especially well suited to spherically symmetric 
space-times. In these coordinates the null geodesic generators of the 
observer’s past light cone are trivially integrated whereas they|-eaBjp.t 
be integrated exactly for Tolman models in (3 -I- coordinates EjUlj. 
The reverse situation holds for the field equations EJ. It follows that, in 
general, the transformation between observational and (3 -I- 1) coordi¬ 
nates cannot be integrated exactly in simple functions. Thus, as could 
be anticipated, the two coordinate systems give different viewpoints on 
the geometry and the physics; the (3 -I- 1) coordinate approach provides 
the physical interpretations and the observational coordinates relate to 
astronomical observations. In this work we needed both viewpoints and 
therefore found it necessary to develop results which enable us to flip 
between the two systems 


2 Notation and Formulae 

In this section a very brief sketch will be given of the leosmological model 
and the notation. For details see Humphreys et al □. In observational 
coordinates the spherically symmetric model takes the forma, 

ds^ = —A{w, y)^dw^ + 2A{w, y)B{w, y)dwdy + C{w, y)^d^ (2) 


where {w = const } label the past light cones on the observer’s world 
line {y = 0 }, j/ is a comoving radial distance coordinate down the null 
geodesic {w = constant, = const } and dQ^ — dO^ -\- s\v? (f)d(f)^. 

The dust velocity is given by u‘ = A~^5l^. The FLRW models are given 
by choosing A = B. In (3 -I- 1) coordinates the model is given by. 


ds^ = —dt^ + 


dR{r,t) ^ dr^ 

dr 1 — kfirY 


+ R{r,tfdn'^ 


( 3 ) 
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where / is arbitrary and relates to the total energy of the system. The 
geometry of the {t = const } hyper-surfaces is parabolic if fc = 0, hyper¬ 
bolic for k = —1 and elliptic for fc = 1. The dust velocity is, u® = SI- 
From the two metrics the subtle relations between the two systems 
can be demonstrated. On a hxed light cone {w = wo} the observer area 
distance (angular diameter distance) D in terms of redshifts is given by, 


= / C{wo,yiz)) 

\R{r{z),T{r{z))) 

where t = T(r) is the equation of the past light cone of the observer. In 
general the function T cannot be determin^ as an exact function and 
so numerical methods have to be employed tj and thus to-determine R 
we need to know both r and t separately. Mustapha et altZI circumvent 
this by restricting attention to a single light cone, which cuts out the 
dynamics, and they use time t which is not an observable. These are not 
appropriate here. 

In the metric (|[) the hyper-surfaces {w = const} are null but that 
is not sufficient to make tlpp the past null cones of the observer. TWe 
requires central conditionstil which develops earlier results of Templets. 
As j/ ^ 0 the following asymptotic behaviour is required. 


A{w,y) = A(w,0) + 0{y) 

( 4 ) 

B{w,y) = B{w,0) + 0(y) 

( 5 ) 

C{w, y) = B{w, G)y + 0{y^) 

(6) 


with A(w,0) yf 0 and B{w,0) yf 0. These conditions imply that near 
enough to the observer the space-time appears Minkowskian. This im¬ 
poses restrictions on the asymptotic behaviour of the number count for¬ 
mula as 1 / rtO. In observational coordinates the number count formula 
is given byliiJ, 

ry 

N{y)=4'K / n{wo,x)B{wo,x)C{'Wo,x)^dx (7) 
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where n(wo, y) is the number density of sources. The dust density is 

p{wo,y) = n{wo,y)m (8) 


where m is the average mass of the sources (galaxies). We ignore selection 
and evolution effects although they could be incorporated without much 
difficulty if we knew their functional form. Also we assume no bias. 
Dark matter could be included via m, if unbiased. From the energy 
conservation equation, the central conditions and the formulae for N 
and p it follows that as y 0, 


P = 


4-kPo 


B{wo,Of 

B(w,0)3 


+ 0{y) 


( 9 ) 
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and 


( 10 ) 


= (I?) 

where po = mn(wo,0) 

These strictly mathematical results mean that the number count for¬ 
mula cannot be fractal for y near zero. 


3 Sketch to Illustrate the Relations Between the Variables 

From the two formulae for the dust velocity we obtain, 


A= — 
dw 

S = -|1 

ay 


( 11 ) 

( 12 ) 


and these can be used to write the held equation for A in the form!) 


A = 


_ C _ 

[(F2 - 1) + (22^)] 


1/2 


(13) 


where C = dC/dw, = FN', F is an arbitrary function of integration 
and a prime denotes the partial derivative with respect to y. We can use 
( ^d] ) to write (§) in the form, 


dt = 


_dC_ 

[(l"^-l) + (2^)]'/' 


(14) 


al ong i t he huid flow held. 

isElS 


Along this how the well known Tolman solution 


dt = 


dR 


[-kf + ^] 


, 1/2 


(15) 


which led us to identify. 


F{yf = 1 - kfiyf (16) 

M{y) = mNt,{y) (17) 

= m [ {1-kf{yf)^^dx (18) 


Clearly for fc = 0, M{y) = mN{y). 

To obtain further equations, e.g. for (1 — kf{y)^), in terms of obser¬ 
vational data D{z) and N{z) requires further use of the held equations. 
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given by Maartens et alu. This results in. 


a 


- kf{zY = 


1+2 


0 ( 2 ) = 1 - 


(1 + x)N'{x) 


D{x)Q{xf 


(l + x)2 


D[x) 


dx. 


dx (19) 

( 20 ) 


The full set of equations is given by Humphreys et altl. The equation 
(@) can be inverted to give, 

1 fy 1 

N{y) = — M’{x)[l-kf{xfr^dx (21) 

^ Jo 


4 An Application of the Formulae 

To apply the formulae we only need to consider one past light cone so we 
can use the area distance D as our distance measure. We will also limit 
ourselves to the k — 0 case which is particularly simple because only 
one of the functions N and D is arbitrary. The other more complicated 
cases which cannsit be handled completely analytically are considered by 
Humphreys et alQ. 

The (impendence of N and D is expressed explicitly by the integral 
equation u, 


1+2 



2mN 

D 


-1 


exp 


m dN 
D'dD 


2mN 

D 


-1 

dD 


( 22 ) 

Once N is known in terms of D, through for instance a fractal formula, 
this gives a relatively simple expression for z in terms of D. Unfor¬ 
tunately the central conditions do not permit such a simple procedure 
since the form of the function N has to change. A simple ansatzu which 
combines the limiting behaviour with a fractal number count formula 
outside the immediate vicinity of D = 0, is 


N{D) 


for D<Di 
D<D<D, 


(23) 


where Dh is the distance at which transition to FLRW geometry occurs. 
We will assume that in the core {D < Di) and in the fractal region 
{Di < D < Dh) the spacetime is parabolic, i.e., k = 0. Then the 
Hubble constant is given bycl. 


Ho = 



(24) 
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which puts a precise constraint on the central density. It can be shownU 
that matching the number count N{y) at Dh relates the four parameters 
Hq, Di, Dh, and v. From observational data one could estimate Hq, Dh 
and v which would then fix the minimum fractal scale Dj. Unfortunately 
the critical observations are controversial. 

For sm^l D an explicit form of the behaviour of 2 with respect to D 
is given byl3, 


dz 

Id 






( 3-^)/2 


vHqDi 




(25) 


which shows that after initial behaviour which is linear by construction 
to Di the z{D) graph curves upwards contradicting the linear Hubble 
Law for scales less than 100 Mpc. This rules out the parabolic models. 

It is interesting that for any rational fractal index v can be 
integrated to give an explicit expression for 2 in terms of D. Note that 
for the simple model of the transition from the core to the fractal region 
given in (M) the integration is particularly simple but it can become 
very complicated for more sophist^ated transition formulae. In the non¬ 
parabolic models it can be shownH that the fractal number count forces 
a very low density for the universe. 


5 Conclusion 

The telk outlined some of the ideas and methods employed by Humphreys 
et alB. Only some aspects have been discussed tonshow the mathemati¬ 
cal ideas involved. The comprehensive treatment □ concludes that exact 
models using fractal counts models out to 150 Mpc either contradict 
the well established linear Hubble law out to 150 Mpc or they yield a 
very low density universe. Both conclusions conflict with observations. 
There are various modifications that could be made to the model but 
the most compelling is that the luminous matter may not trace the ac¬ 
tual distrilmtion of density. In this connection a recent paper by Labini 
and DurertJ argues for a fractal baryonic matter distribution with non 
fractal dark matter. The approach developed here gives an alternative 
method by which bias could be detected or measured. 
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